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1. Introduction 

Estimates from the current pulsar birth rate and 
the number of supernovae to account for the 
heavy-element abundance suggest that in total 
10 8 to 10 9 isolated neutron stars exist in our 
Galaxy. Only a small fraction is detectable as 
young neutron stars due to their emission of 
thermal X-rays (for ~10 6 years) or pulsed ra- 
dio emission (for ~10 8 years). Proposals that 
an old "recycled" population of isolated neu- 
tron stars re-heated by accretion from the in- 
terstellar medium (ISM) should be detectable 
by ROSAT, triggered several projects to search 
for such objects in the ROSAT all-sky sur- 
vey data. Over the last decade seven very soft 
X-ray sources with particular characteristics 
were discovered in ROSAT data. Extreme X-ray 
to optical flux ratios and low absorption col- 
umn densities strongly suggest that these ob- 
jects are nearby isolated neutron stars (see re- 
views by iTreves et al.ll2000t lMotchll2001t lHaberll 
I2nn4fi . Using HST, parallax measurements yield 
a distance of 117±12 pc for RXJ1856. 5-3754 
{Walter fc LattimerfEio^ . The detection of rel- 
atively high proper motion (PM) for the three 
brightest stars makes accretion from the ISM 
highly ineffective and favours the picture of cool- 
ing neutron stars with an age of ~10 5 -10 6 years 
to power the X-rays. Tracing back the apparent 
trajectories suggests that the brightest of the 
ROSAT-discovered isolated neutron stars were 
born in the Sco OB2 comple x which is the clos- 
est OB association (see e.g. iMotch et al.ll2005l 
and references therein). 

The X-ray spectra of the "magnificent 
seven", as they are sometimes called in the lit- 
erature, are thermal and blackbody-like without 
a hard power-law tail as it is ofte n observed in 
other isolated neutron stars (e.g. iPavlov et alJ 
2002). Typical observed blackbody tempera- 



tures kT are in the range of 40—110 eV (see 
Table 1). From five stars X-ray pulsations were 
detected with pulse periods between 3 s and 
12 s and pulsed fractions between 4% and 18% 
(Fig.Q. However, for the X-ray brightest star 
RX J1856.5— 3754 no pulsations with a strin- 
gent upper limit on periodic variation of 1.3% 
(2cr confiden ce level in the 0.02 - 1000 s range) 
were found {Burwitz et alJ 12003ft . A surprising 
discovery was that the X-ray spectrum and the 
pulsed fraction observed from RX J0720.4— 3125 
changes on a time-scale of years which may 
be caused by prece ssion of the neutron star 
(|de Vries et al.ll2004ft . 

2. Broad absorption lines 

XMM-Newton observations of the thermal 
isolated neutron stars revealed deviations from 
the Planckian shape in the X-ray spectra 
obtained by the EPIC-pn and RGS instru- 
ments. Fig.|3 shows a comparison of the 
EPIC-pn spectra of the six best observed 
thermal isolated neutron stars fitted with 
an absorbed blackbody model. Large resid- 
uals are seen from RB S 1223 (iHabe rl et al.l 
12003ft . RXJ0720.4-3125 {Haberl et alJ l2004bh 
and RXJ1605. 3+3249, in the latter case 
the devi ations were discov e red in RGS 
spectra {van Kerkwiik et all 120041) . Non- 
magnetic neutron star atmosphere models (e.g. 
lOansicke et all2002HZavlin fc Pavlml2002ft can 
not explain the X-ray spectra: Iron and solar 
mixture atmospheres cause too many absorp- 
tion features and deviations from a blackbody 
model in particular at energies between 0.5 and 
1.0 keV which are not seen in the measured 
spectra. On the other hand the spectrum of 
a pure hydrogen model is similar in shape to 
that of a blackbody and does not fit the data 
either. Moreover, hydrogen atmosphere models 
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Table 1. X-ray and optical properties of the magnificent seven 



Object 



kT Period Amplitude Optical PM 
eV s % mag mas/year 



Ref. 



RXJ0420.0-5022 


44 


3.45 


13 


B = 26.6 




1 


RXJ0720.4-3125 


85-95 


8.39 


8-15 


B = 26.6 


97 


2,3,4,5,6 


RXJ0806.4-4123 


96 


11.37 


6 


B > 24 




7,1 


RBS 1223 (1) 


86 


10.31 


18 


ni50ccd = 28.6 




8,9,10,11 


RXJ1605.3+3249 


96 






B = 27.2 


145 


12,13,14,15 


RXJ1856.5-3754 


60 




<1.3 


V = 25.7 


332 


16,17,18 


RBS 1774 (2) 


101 


9.44 


4 


R > 23 




19,20 



(1) = 1RXS J130848.6+212708 

(2) = 1RXS J21 4303. 7+06541 9 

Referen ces: (11 lHaberl et all (l2004aD (2) lHaberl et alJ (119971) (3) ICropper et al 



2004 

1999]. 

19991) (13) iKaplanetafl 



fcOOll) ( 4) lHaberl et al 



(5 ) Ide Vries et al.1 l20oiTO d) iMotch et all J2003D (7) lHaberl fc Zavlinl p002\) (8) ISchwope et al 
f 9) lHambarvan et alJ <2002ft flO) iKaolan et alJ i2002ft fi ll lHaberl et alJ <200fl fl2TlMotch et al 



.. ._ (14) Ivan Kerkwiiket all (l2004fl i (15) lMotch et al.1 120051) (16) 

alter fc Matthews! <199ffl (1 7) rWalter fc Lattimerl (^2002f) (18) iBurwitz etafl lj2003l) f 19) IZampieri et ail 



(20) IZane et alJ l|2005D 



Table 2. Magnetic field estimates 



Object 


dP/dt 


E C y C 


Bdb 


B C y C 




lO" 13 ss" 1 


eV 


10 13 G 


10 13 G 


RXJ0420.0-5022 


<92 


330? 


<18 


6.6? 


RXJ0720.4-3125 


1.4±0.6 


260 


2.8-4.2 


5.2 


RXJ0806.4-4123 


<18 




<14 




RBS 1223 


<9 


100-300 


<10 


2-6 


RXJ1605.3+3249 




450-480 




9.1-9.7 


RXJ1856.5-3754 






~1 




RBS 1774 




-700 




~14 



over-predict the actually obser ved optical fluxes 
by l arge factors (^300, see IZavlin fc Pavlov! 

l2onl . 

The XMM-Newton spectra can best be mod- 
eled with a Planck continuum including a broad, 
Gaussian shaped absorption line (Fig.[3J). Line 
centroid energies are summarized in Table 2. 
In the EPIC-pn da ta of RBS 1223 (see also 
ISchwope et alJl200Hl) and RX J0720.4-3125 the 
depth of the absorption line (or the equiv- 
alent width) was found to vary with pulse 
phase. In the cases of RX J0806.4— 4123 and 
RXJ0420.0-5022 it is not clear to which 
extent the residuals of the blackbody fits 
are caus ed by systematic ca libration uncer- 
tainties l|Haberl et alJ l2004a|b In particular 
RX J0806.4— 4123 shows a residual pattern sim- 
ilar to that of RXJ1856.5-3754 which is be- 
lieved to exhibit a pure blackbody spectrum as 
seen from the high resolution C handra LETGS 
spectrum ijBurwitz et al.ll2003[) . For RBS 1774 
recently an absorption feature at ~0.7 keV was 



reported from th e analysis of EPIC spectra 
l|Zane et al.ll200filb At such high energies (the 
highest line energy reported from the thermal 
isolated neutron stars) the energy resolution is 
better and the calibration uncertainties smaller 
than at lower energies. 

3. Strongly magnetized neutron stars 

An H Q emission line nebula was dis- 
covered around RX J1856.5— 3754 

l)van Kerkwiik fc KulkarnH l2001j) . Assuming 
that magnetic dipole breaking powers this neb- 
ula and using an age of th e star of 5 x 10 5 years 
l|Walter fc Lattimerll2002T) allows an estimate of 
the magnetic fi eld strength of the neutron star 
of B ~ 10 13 G llBraie fc Roma,nill2nn2HTriimnerl 
2005). A similar magnetic field strength of 
(2.8—4.2) x 10 13 G was inferred from the pulse 
period history of RX J0720.4— 3125 as observed 
with ROSAT, Chandra a nd XMM-Newton ove r 
a time span of 10 years l)Cropper et al J 12004]) . 
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EPlC-pn (0.12-1.2 keV) Rev 0534 RX J0720. 4-3125 EPIC-pn (0.12-1.2 keV) Rev 0168 RX J0806.4-4123 





0.5 1 

Phase 

EPIC-pn (0.12-1.2 keV) Rev 0377 



1RXS J130848. 6 + 212708 



0.5 1 

Phase 

EPIC-pn (0.12-0.7 keV) 



1.5 

RXJ0420. 0-5022 





Fig. 1. EPIC-pn light curves folded with the pulse period of four thermal isolated neutron stars. 
For direct comparison of the pulse profile the flux is normalized to the mean and plotted on the 
same scale. Except for RX J0420. 0—5022, which shows the softest X-ray spectrum, data from the 
same en ergy band was used. To gain statistics data of four observations of RXJ0420. 0—5022 were 
merged l|Haberl et al.ll2004al) . 



These were the first indications that the group 
of thermal isolated neutron stars possess strong 
magnetic fields of the order of 10 13 — 10 14 G. 
Such strong fields are indeed required to spin 
down the neutron stars to their current long 
rotation periods within 10 5 — 10 6 years (still 
being sufficiently hot to be detected in X-rays) 
if they were born with msec periods. 

Cyclotron resonance absorption features in 
the 0.1 — 1.0 keV band are expected in X-ray 
spectra from magnetized neutron stars with field 
strengths in the range of 10 10 — 10 11 G or 2 
x 10 13 — 2 x 10 14 G if cause d by electrons or 
protons, respectively ( see e.g. IZane et al"ll200ll 
IZavlin fc Pavlov| 1200 2'). Variation of the mag- 
netic field strength over the neutron star surface 
(as expected for d ipole fields) lead s to a broad- 
ening of the line l)Ho fc Laill2004[) . The strong 
magnetic fields inferred from magnetic dipole 



breaking effects in at least two of the stars sug- 
gests that the broad absorption features seen in 
the X-ray spectra of thermal isolated neutron 
stars originate from cyclotron resonance absorp- 
tion by protons or highly ionized atoms of heavy 
elements. With a mass to charge ratio of ^2 with 
respect to protons the latter case would lead to 
B a factor of ^2 higher than that derived for 
protons. Different ionization states would result 
in a series of lines with energies differing by only 
a few percent, leading to additional broadening 
of the lines. 

An alternative possibility for the origin of the 
absorption line is atomic bound-bound transi- 
tions. In strong magnetic fields atomic orbitals 
are distorted into a cylindrical shape and the 
electron energy levels are similar to Landau 
states, with binding energies of atoms signifi- 
cantly increased. E.g. for hydrogen in a mag- 
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Fig. 2. EPIC-pn and RGS spectra of thermal isolated neutron stars fitted with an absorbed black- 
body model. The fits represent the calibration status as available with SAS release 6.0.0. 



netic field of the order of 10 13 G the strongest 
atomic transition is expected at energy E/eV 
75 a+0.131n(Bi 3 ))+63Bi3 , with B 13 = B/10 13 
G jZavlin fc Pavlovll2002ft . For the line energies 
found in the spectra of thermal isolated neutron 
stars this would require similar field strengths 
to those derived assuming cyclotron absorp- 
tion. Atomic line transitions are expected to be 
less prominent at higher temperatures because 



of a higher ionization degree l(Zavlin fc Pavlov! 
120021) . 

4. Conclusions 

Although the true origin of the broad absorption 
lines in X-ray spectra of thermal isolated neu- 
tron stars is not clear yet, our current knowledge 
about the "magnificent seven" strongly suggests 
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Fig. 3. Top two panels: EPIC-pn and RGS spec- 
tra of two thermal isolated neutron stars fitted 
with an absorbed blackbody model including a 
broad absorption line. Bottom panel: RGS spec- 
trum of RX J1605.3+3249 modeled with black- 
body continuum and broad absorption line. A 
possible additional narrow absorption line at 
21. 5A is also included ffrom lvan~K crkwii k et al.l 
12004) . 



that they are highly magnetized (10 13 — 10 14 G), 
slowly rotating cooling neutron stars. Further 
timing studies would be very useful to obtain 
more independent estimates of the magnetic 
field strength (as they currently only exist from 
RXJ0720.4-3125). 

We do not detect radio emission, prob- 
ably because their radio beam is very nar- 
row due to their large light cylinder radius. 
The discovery of a few radio pulsars with 
simila r magnetic fi e ld str e ngth and long pe- 
riod iCamilo et alJ 120001 iMorris et aD 120021 
lMcLaughTm*^rii,lJl2003|) shows that radio emis- 
sion can still occur at inferred field strengths 
higher than the "quantum critical field" B cr — 
m 2 e <?/eh ~ 4.4 x 10 13 G. On the other hand, 
any non-thermal emission from the "magnificent 
seven" may so far just fall below th e detection 
thres hold of current instruments ijZane et alJ 
l2005|) . 
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